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Introduction 


The numerical treatment of the Schrédinger equation in quantum chemistry 
gives rise to molecular integrals of different kinds. The solving of these integrals 
is one of the most tedious problems in determining the energy levels of a certain 
molecule. Especially for crystals, where a large amount of atom neighbours must 
be taken into account, it is necessary to have a rapid way of computing these 
integrals. 

For this purpose some of the most common one-electron integrals have been 
programmed by the present author for the Swedish electronic computer Besk. 
This programme is completely automatic and contains one-electron two-center inte- 
grals for all combinations of Slater atomic orbitals of the K, L and M shells. For 
integrals containing orbitals of higher electronic shells the programme is not appli- 
cable. Many of these integrals can, however, be rather easily computed by using 
the molecular functions A,(«) and B,(6), which have been tabulated by the present 
author for n<8[1]. They can be given for a continuous range of values of the 
parameters « and f by the present programme and they are useful also for the 
computation of two-electron integrals. In an earlier note [2], the author has pre- 
sented a programme for overlap integrals. This programme is included as a special 
case in the present programme. 

Many excellent tables of molecular integrals have already been made by several 
other authors (e.g. Kopineck [3], Barnett and Coulson [4], Kotani et al. [5], Roo- 
thaan [6], Preuss [7], Miller, Gerhauser and Matsen [8] and many others). These 
tables also give other types of molecular integrals than those of the present pro- 
gramme. The Chicago group (Roothaan et al.) has also made extensive programmes 
for electronic computers containing all types of one-electron integrals and two-elec- 
tron one-, two- and three-center integrals. All these integrals, however, contain 
only K and L shell orbitals. Thus the inclusion of M shell orbitals in the present 
programme will give some new material in addition to other available tables and 
programmes for electronic computers. The present programme has now been tested 
and used frequently for more than two years. Thus the probability that there still 
are cases for which the programme might fail is rather small. 
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Fig. 1. Co-ordinate systems for the two nuclear centers a and b. 


1. Description of the programme 


The following types of one-electron integrals with two centers, a and b, can be 
computed by use of the programme: 


A: Overlap integrals, | yr WL ite 


. a 1 b a 1 b 
B: Resonance integrals, [ WE Prat; | pr— pr dt. 
J ra J rp 
ae é : a it b 
Kinetic energy integrals, [ Wr (—ZA) yr dt. 


a 


. . a 1 a 
D: Nuclear attraction integrals, [ Wr—Yprdt. 


rp 
The molecular A and B functions 
ce 1 
A, (x)= \ ee dé, B,(B)=| ne PI dy 
i 
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can also be computed by use of the programme. 

The orbital indices AK and L! are used for indication of the integrals, and represent 
different nodeless hydrogen-like atomic orbitals of Slater type, according to the 
following table: 


ig (0) 1 2 3 4 5 6 a 8 is) 
we ls 28 2p% 2G se 8p Boo 38do Sdm 3do 


The nuclei are denoted by a and b, and the co-ordinate systems for the corre- 
sponding two centers have been chosen according to Fig. 1. We note that the two 
o-directions (z, and z,) are chosen to be parallel. This means, for example, that the 
overlap integrals (13) and (31) have opposite signs. The same thing is true for 
resonance and kinetic energy integrals. 


. 


2. Use of the programme 


The present programme is very simple to use and is available for everybody, 
who is interested in using it for computing molecular integrals. 

There are three numerical parameters, which must be given for each integral, 
namely (cf. [2]) 


1 These indices do not correspond to the K and L shells. 
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R =half the inter-atomic distance 
Zx = effective charge of the K orbital 
Z, = effective charge of the L orbital 


All these data values must be given in atomic units and put on the parameter tape 
in decimal form in the sequence given above and finally followed by the notation A. 
The decimal data values are written without decimal point and every value is fol- 
lowed by its sign notation (D for plus, B for minus) and the 10-exponent of its 
scale factor also with its sign notation (see example below). 

After the data values the pertinent notation KL for the participating atomic 
orbitals and the notations ABCD for the different types of integrals are written. 
If some of the integrals are not wanted, it is possible to exclude some of these 
notations. However, since the kinetic energy integrals, C, include linear combina- 
tions of A and B, the latter notations must always precede the notation C and 
cannot be excluded in order to get C correctly computed. For similar reasons, none 
of the notations A, B or C can be written after D. 

The notations ABCD are finally followed by the notation #, and another para- 
meter scheme can be added. In the last parameter scheme, H must be replaced by 
fF, which causes Besk to stop. 

If the molecular A and B functions are wanted, the parameter scheme can be 
written as follows. At first the numerical data values of « and f, resp., must be 
written according to the decimal data rule, followed by the code notation A. Fi- 
nally the notation AB is written. Then the scheme is finished and another one can 
be added. There is no certain final stop notation for this scheme. 

The tape of the Besk programme, which has the notation S83, is available at 
Matematikmaskinnimnden in Stockholm. Together with the parameter tape, writ- 
ten according to the scheme above, it causes Besk to print an output which is de- 
scribed below. 

The actual notation K LZ is printed to the left in the output and is followed by 
the data values R, Z; and Z,, written as ordinary decimal fractions. The values 
of the different integrals then follow, preceded by their pertinent notations. They 
are given with eight decimals together with the 10-exponents of their scale factors 
(see below). The notation B is followed by both the types of resonance integrals 
in the sequence given above. 

As an example, we suppose that the B and D integrals of (28) are wanted for 
R=1.6251, Z,=2.41 and Z,=1 and the same integrals of (16) for k=3.9233, 
Z, = 2.68 and Z,=1.3 (a.u.). The parameter tape must then be stamped as follows: 


1G251D1ID: 2P4AID1D -1DID A 28 © BD OE 
2003301) 26801) 13DiID A 16 BD PF 


The corresponding output will be 


28 1.6251 2.41 1 B —07025769 00 —04273470 00 D 09261602 -01 
16 3.9233 2.68 1.3 B — 14653864 00 —0416152700 D 06612444 -01 


This means that the first resonance integral of (28) (i.e. for K=2, L=8) for the 
given parameters is equal to —0.07025769 a.u., etc. 
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Fig. 2. Diagrams of overlap integrals, A, for orbitals with equal effective charges, Z (= principal 

duaneum number times orbital exponent in a.u.), as functions of ZR (R=half the inter-atomic 

distance a—b in a.u.). The different curves are labelled by the two involved orbital indices 
K and L. 


In the output for the molecular A and B functions the notation AB is printed, 
followed by the values of « and #. The A and B functions are then given in two 
tables, the left one a A,(a) pnd the right one B,(8) for n=0,1,... 8, if |p | 22, 
but only for n=0, 1,... 6, if |B|< 


3. Diagrams for some integrals of the K, L and M shells for equal exponents 


As mentioned in section 2, the three parameters R, Z, and Z, are needed for the 
description of the two-center integrals. However, the overlap integrals, A, will 
remain unchanged, if R is divided by an arbitrary constant, k, at the same time 
as Z, and Zy are multiplied by k. In tabulating these integrals, thus only two 


parameters are needed. For example, these parameters can he chosen as RZ, and 
RZ,, or as one of these values and the ratio Zy/ Za: 
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T 


Fig. 3. Diagrams of B/Z as functions of ZR (equal Z). B =resonance integral (a.u.). Integral 
curves of the K, L and M shells have been plotted. 


The other types of integrals are not invariant under transformations of this kind. 
However, there are certain simple connections between integrals transformed in this 
way, which makes it unnecessary to use more than two parameters for tabulating 
them. Choosing the constant k equal to FR, we can write these connections for B 


and D integrals as 


1 
B (BR, Zq, Zy) = RP (Lele Za ly 5) 


and for C integrals as 


1 
CO (R, Za, Zo) = Fa OC 1, RZ, B Z,). 


The diagrams given here correspond to equal effective charges, 1.e. Z,= Z) = Z. 
Consequently, mainly diagonal integrals have been given. These integrals can be 
described by diagrams with only one parameter, RZ. The integrals for Z=1 are 


found directly from the diagrams. The remaining integrals with equal effective 
charges, can be found by using the following formulae 
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The Figs. 2-5 show different diagonal integrals of the A, B, C and D types (for 
Z=1) as functions of the variable RZ. To find the corresponding integrals for 
Z=1, the B and D values for Z=1 must accordingly be multiplied by Z and the 
CO values by Z, while all the A integrals are found directly. Some non-diagonal over- 
lap integrals combining equal and neighbouring shells have also been given, although 
these integrals might rather have been given with different effective charges. 

It is interesting to notice the well-known fact that the different shells roughly 
correspond to different characteristic interaction ranges, which, however, also de- 


86 


ARKIV FOR Fysik. Bd 17 nr 3 


Oo 
sete) 
N 


oO 
h 
al 


0.20 


+44 + 4 +4 +4 ++ 


0.05+ 
+ 
+ 

0 — ~ —+- + a ~ + . +— + + + ZR 


Fig. 5. Diagrams of D/Z as functions of ZR (equal Z). D =nuclear attraction integral (a.u.). 


pend on the orbital exponents. By suitable choices of the orbital exponents these 
ranges can be made roughly equal for different shells. For example, if the ratio 
Zu/Z, of the Z-values for the M and L shells, resp., is chosen equal to about 2, 
this is seen to be the case for the L and M shells from Fig. 2. For physical reasons, 
however, the effective charges of the M shell orbitals of an atom would be less than 
those of the LZ shell orbitals, since M electrons would be more screened by the inner 
electrons than JL electrons. Nevertheless, it might be of value to consider also 
the opposite case since physical pictures based on the conditions in the free atom 
are not always useful. These pictures may be applicable to atoms and small mole- 
cules, but for crystals they must not always be taken too literally. The use of 
Slater exponents for atomic orbitals of excited shells usually gives very large over- 
lap for the nearest neighbours in a crystal and, consequently, bad convergence for 
the energy terms with increasing number of neighbours. Thus, the same orbital 
exponents as those found for atoms and molecules cannot always be used also for 
crystals. Actually, variation shows that much larger values often have to be used. 

The diagrams given here are useful for finding the upper limit of atomic neigh- 
bours, which ought to be taken into account in energy band calculations for 
crystals. 

As is seen from the diagrams, the shapes of the potential and kinetic energy 
integrals are rather similar to the shapes of the corresponding overlap integrals. 
However, they are far from proportional and their ratios sometimes even change 
their signs. Especially the kinetic energy integrals, C’, differ rather much from being 
proportional to the overlap integrals, since they decrease more rapidly (for small Z 
values) with increasing argument, ZR. The functions D/Z all decrease as the in- 
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verse inter-atomic distance, 1/(2R), as is seen from Fig. 5. Consequently the ap- 
proximation of considering overlap and energy elements as being proportional, 
frequently used in quantum chemistry, is far from satisfactory. However, using a 
suitable effective Hamiltonian operator for the system under consideration, the 
one-electron energy integrals given by the present programme, are useful for a rough 
calculation of the electronic states. 


8. 
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